The stability of Killing-Cauchy horizons 
in colliding plane wave space-times 



in 
o 
o 

C 

a 
in 



o 



J. B. Griffiths* 

Department of Mathematical Sciences, Loughborough University, 
Loughborough, Leics. LE11 3TU, U.K. 



February 7, 2008 



Abstract 



It is confirmed rigorously that the Killing-Cauchy horizons, which sometimes occur in space- 
times representing the collision and subsequent interaction of plane gravitational waves in 
a Minkowski background, are unstable with respect to bounded perturbations of the initial 
waves, at least for the case in which the initial waves have constant aligned polarizations. 



O ■ 1 Introduction 

in ' 

Many classes of explicit exact solutions are known which model the collision and subsequent 
interaction between shock-fronted plane gravitational waves which propagate and collide in a 
i ! Minkowski background (for reviews see CP or [2]). In all these solutions, some kind of singularity 
£jpj! always appears in the interaction region. This is generally a spacelike curvature singularity (like 
the time-reverse of an initial cosmological singularity). However large classes of solutions with 
an infinite number of parameters exist in which the scalar polynomial curvature singularity is 
replaced by a horizon. In these cases, the space-time can be extended through the horizon, but 
the extension is not unique. It has become widely believed that such Killing-Cauchy horizons 
are unstable with respect to small changes in the initial data. For example, Yurtsever |3j has 
shown that they are unstable with respect to the addition of some perturbative linear field which 
preserves the G2 symmetry. However, this result does not answer the question of whether they 
are unstable with respect to variations in the approaching purely gravitational waves when the 
vacuum field equations are satisfied exactly. It is the purpose of the present paper to investigate 
this question in detail. 

It must first be pointed out that, in the vast amount of work that was undertaken on 
this topic in the 1970s and 1980s, and in all the exact solutions that were then produced, the 
approach was adopted of first solving the field equations in the interaction region. Once a family 
of such solutions had been obtained, the free parameters were constrained to satisfy the junction 
conditions that have to be imposed in order to extend the solution to the prior regions. Thus, 
the approaching waves which physically give rise to the solutions were only determined once 
the solution had been obtained. Within this context, it was argued that those solutions which 
contained horizons were unstable with respect to perturbations of the initial data. However, it 
was not then appreciated that the perturbations which transform a Killing-Cauchy horizon to a 
scalar polynomial curvature singularity also normally introduce singularities in the initial waves 
prior to their collision. The question therefore still needs to be addressed as to whether or not 
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the horizons are stable with respect to regular (bounded) perturbations of physically acceptable 
initial waves. 

It is only very recently that techniques have been developed by which colliding plane wave 
solutions can be explicitly constructed from their characteristic initial data representing the 
approaching waves [I]-[E]- It is therefore only now that tools are available to reconsider this 
question. In this paper, only the linear vacuum case in which two approaching gravitational 
waves have constant aligned polarizations will be considered. This situation is much easier to 
analyse, but it already demonstrates the essential features of the physical situation. The basic 
result, not surprisingly, is that the Killing-Cauchy horizons that sometimes appear in colliding 
plane wave space-times are unstable with respect to bounded perturbations of the initial waves 
that generate these solutions. 



"Focussing" singularity 




Figure 1: The structure of colliding plane wave space-times. Region I is the background Minkowski 
space, regions II and III contain the approaching plane waves, and region IV represents the interaction 
region following the collision. The singularity structure following the collision is known from the study 
of families of exact solutions and is described in £Q. 



2 Initial data 

Colliding plane wave space-times are naturally divided into four regions as indicated in figure 1. 
It is found to be convenient to use two future-pointing null coordinates u and v throughout the 
space-time. The four distinct regions can then be identified as those that are separated by two 
null hypersurfaces taken as u = and v = which represent the wavefronts of the approaching 
waves. 

The background region I (u < 0, v < 0) is a flat vacuum represented by the line element 

ds 2 = 2 du dv - dx 2 - dy 2 . (I) 

Region II (u > 0, v < 0) contains one of the approaching plane waves. If this wave has 
constant (linear) polarization, it can be described by a metric in the Brinkmann form 

ds 2 = 2 dudr - dX 2 - dY 2 + h + {u){X 2 - Y 2 ) du 2 , (2) 

where h + (u) = x I , 4 + (n), which is the only non-zero component of the Weyl tensor (relative to 
an appropriate tetrad). This arbitrary function explicitly represents the profile of the plane 
gravitational wave as a function of the retarded time u. (The subscripts + and — are used here 
and below to denote functions that are defined in regions II and III respectively.) 

The form of the metric © may be taken to include region I with h + = for u < 0. However, 
to analyse the collision of plane waves, it is appropriate to initially transform the line element 



PJ) to the form 

ds 2 = 2dudv- e 2 ^+ dx 2 - a\ e~ 2 ^+ dy 2 . (3) 

As explained in jSj, the coordinate it is an affine parameter on the null geodesies v = const. 
x = const, y = const, which cross the wavefront from the background region. Also, the 
function ip+(u) must satisfy the linear scattering equation 

(e*+\ m + h + (u)e*+=0. (4) 

In order to join (J3J) smoothly with (J2), it is necessary that ^+(0) = and that ip +tU (0) must 
vanish except for any non-zero component which arises from a possible impulsive component on 
the wavefront. With these conditions, ip+(u) is determined uniquely from the initial wave profile 
h + (u). The function ip + (u) can therefore be taken to represent the initial data to be prescribed 
on the characteristic v = 0. 

With ip+(u) determined, the initial gravitational wave is given by the Weyl tensor component 

^4+(u) = -tp +)U S(u) - [ip+, uu + (V>+,«) 2 ] <3>(u). (5) 
The other metric function a + (u) in Q is then determined for region II from the linear equation 

a+, uu ~ 2 ip+, u a+, u + 2(i/j +tU ) 2 a+ = 0, (6) 

whose solution is uniquely specified by the initial data a+(0) = 1, and a +jU (0) = 0. It can be 
seen that a + (u) must be a monotonically decreasing function in this region. 

Region III (u < 0, v > 0) contains the other plane wave which approaches from the opposite 
direction. If this wave also has constant polarization which is aligned with that in region II, the 
metric can simply be taken as having the same form as that in that region, but with the roles 
of u and v reversed. Thus, the initial data that is specified in this region will be taken as the 
function ip_(v). The non-zero component of the Weyl tensor *&q_(v) representing the second 
initial gravitational wave is given by an equation equivalent to (JSJ), and the metric function 
a^(v) is given by the equivalent of © and the corresponding initial data. 

Region IV (u > 0, v > 0) represents the region in which the waves interact. In the vacuum 
case for the collision of plane gravitational waves with constant aligned polarization, this region 
can always be described by the line element 

ds 2 = 2fdudv-e 2 ' 4 'dx 2 -a 2 e- 2i 'dy 2 , (7) 

where ip(u,v), a(u,v) and f(u,v) are now functions of both null coordinates. If characteristic 
initial data is taken as described above, these functions must satisfy the initial data: 

ip(u,0) = ip+(u) a{u,0) = a+(u) f(u,0) = 1 , . 

ip{0,v) =tp_(v) a(0,v) = a_(«) f(0,v) = l U 

with ^(0,0) = 0, a(0, 0) = 1 and /(0, 0) = 1. The necessary (vacuum) field equations will be 
considered in the following section. 



3 The field equations 

One of the vacuum field equations for the metric is 

whose solution, satisfying the necessary initial conditions, is given by 



a(u, v) = a + (u) + a + (u) — 1. 



0) 



It is convenient to introduce the functions £(u) and i]{v) such that 2a = £(u) — rj{v), where 



, , j 1 in regions I and III . , J — 1 in regions I and II 

I 2a + {u) — 1 in regions II and IV ' ' ' j 1 — 2a_(i>) in regions III and IV 

Then, since a + (u) and a_(i>) are monotonically decreasing, £(u) must be a monotonically de- 
creasing function and rj(v) monotonically increasing in the interaction region. 
The other main vacuum field equation is 

2 a ip iUV + a tU ip >v + a „ ip >u = 0, (10) 

where a is given by Q and the initial data by ©• For any solution of IjlOJI satisfying the 
initial data, the additional metric function / can be obtained by integrating the remaining field 
equations which are 

f )V a vv 2a^ 2 /„ a ;UU 2aip tU 2 

-j- = — — i ! — 2y,v, -j- = — — i ! — 2y>,w (if) 

The non-zero components of the Weyl tensor in the interaction region are then 

^0 = -i>,vO{v) ' j— 1 1 j^- 

^2 = (2a ip, v i> )U - a jU ip tV - a >v ip, u ) (12) 

^4 = -1p,u0{u) ! ! h 



/ p 

It may be observed that a solution can be determined in region IV up to the spacelike hypersur- 
face on which a = — rj) = 0, on which some kind of singularity will occur. This is referred 
to as a "focussing" singularity in figure 1. Its character will be discussed below. 

In view of the properties described above, it is always possible to adopt £ and rj as coordinates 
throughout the interaction region, in which £ < 1, r/ > — 1 and £ — 77 > 0. In terms of these 
coordinates, the main field equation (jTU|) becomes 

(Z-V)^- 5^ + 5^ = 0, (13) 

which is an Euler-Poisson-Darboux equation with non-integer coefficients, whose solution rf) 
must satisfy the initial data 

^,-i) = M0, ^0-,v) = Mv), (14) 

with —1) = 0. This is sufficient to uniquely determine the solution of the colliding plane 
wave problem in the interaction region. However, as explained in [Sj, to satisfy the initial 
data, if}(£, rf) must be nonanalytic on each wavefront where it behaves like yT — £ or y/1 + r] 
respectively. Also, in terms of these coordinates, the non-zero Weyl tensor components are 

= ~V,v tv s ( v ) - % (tm + 2(£ " V) tv 3 + W/) 



*-> - (2(£ - v)i>* +tt - tv) ( l5 ) 

S.u. 
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*4 = ~Z,u t^(u) - - 2(e - V ) 4>/ + W/) • 



4 The singularity in the interaction region 



Now consider the singularity which occurs in region IV on the spacelike hypersurface on which 
a = l(£ — 77) = 0. It can be shown (see references in pQj ) that this is generically a scalar 
polynomial curvature singularity on which the invariant ^o^4 — 4 x I'i^'3+3^ / 2 diverges. However, 
families of solutions for colliding plane waves exist, even having infinitely many parameters, in 
which this invariant does not diverge. In these cases, the singularity at a = corresponds to a 
Killing-Cauchy horizon through which the space-time can be extended. 

The general solution of the main equation (|13|) can be expressed (exactly as for vacuum 
Gowdy cosmologies) in terms of the coordinates a = — 77) and (3 = + 77) in the form 

e^duj, (16) 

where Jq{ujq) and Yq[lool) are Bessel functions of the first and second kinds of zero order, and 
a(u) and b(u) are arbitrary functions of the parameter u>. The terms which involve Jq{ujo) 
remain well behaved on the singularity at a = 0, while the terms which involve Yq(uo) diverge 
logarithmically (the presence of some terms of this type is necessary to satisfy the initial data). 

In his analysis of the asymptotic behaviour of colliding plane wave space-times, Yurtsever 
[S] has shown that, for any timelike geodesic which approaches the singularity a = 0, the metric 
asymptotically approaches that of a particular Kasner solution. However, the particular Kasner 
exponents vary for geodesies which approach different points on the singularity. It follows that 
the focussing singularity is generally a curvature singularity. The only exception is that in 
which the metric approaches that of the degenerate (flat) Kasner solution at all points over 
the hypersurface a = 0. In this case, the focussing hypersurface is a Killing-Cauchy horizon 
across which space-time can be extended. Such solutions are characterised by the asymptotic 
behaviour 

7/^,77) -log a -£(/?) or 7/^,77) ~B(f3) (17) 

as a — ► 0, where B(f3) is an arbitrary function which is independent of a. However, the 
substitution ip — > log a — ip simply corresponds to an interchange of the x, y coordinates. Thus, 
the two possible conditions in Q17JI are physically equivalent, and the second can be adopted as 
the criteria for a colliding plane wave solution with a horizon. 

In considering perturbations of such solutions, Yurtsever (§| has rigorously shown that ar- 
bitrarily small changes in the amplitude functions (specifically in the coefficients b(u>) in (|16j) ) 
will cause variations in the Kasner parameters. He has reasonably concluded that, although 
there are colliding plane wave space-times which contain a Killing-Cauchy horizon rather than 
a space-like curvature singularity, these space-times are unstable against small perturbations of 
the initial data and that 'generic' initial data leads to curvature singularities. However, it may 
be observed that the perturbations which induce a curvature singularity in the interaction region 
also generally have the effect that the relevant Weyl tensor component becomes unbounded at 
the fold singularities in the initial regions II and III. In other words, the components which 
introduce a curvature singularity into a space-time which otherwise has a horizon are associated 
with unphysical initial data. Thus, the question that still needs to be addressed is that of the 
stability of the horizon with respect to perturbations which correspond to bounded variations 
of the initial gravitational waves prior to their collision. 



tp = a(u))J (uja) + b(uj)Y 



5 Explicit solutions with horizons 

Now consider the colliding plane wave solutions jlUj for which 

* = log a - ^ cosh- 1 - k - cosh- 1 {l±£) , 



where a, k + and k_ are constants. When a = and k + = k_ = 1, this is the Khan- 
Penrose solution jllj which represents the collision of purely impulsive waves. It also includes 
the degenerate Ferrari-Ibahez solutions for which the interaction region is locally isomorphic 
to part of the Schwarzschild solution inside the horizon. 

By showing that, as a — > 0, the scalar invariant behaves as 



*o^4 + 



2 8 / 2 



■((a + k + + k_) 2 - l) 2 ((a + k + + k„f + 3)a- 



Feinstein and Ibahez {T3] have shown that colliding plane wave solutions given by (|18|) have a 
horizon at a = in region IV rather than a curvature singularity if a + k + + k^ = dbl. (This 
also follows from the condition (|17)l.) 

It can be seen that the last two terms in the expansion (|18|) are required to describe the 
necessary behaviour of the solution near the wavefronts £ = 1 and rj = — 1 respectively, on which 
they also become zero. Also, all three components of (|18[1 diverge logarithmically as a — > 0. 
Thus, these terms must include components which arise from the Yq{u)o) terms in (|16|) together 
with possible components which arise from the regular terms. In addition, it is clear that the 
solutions (|18jl arise from the initial data functions 



' i 



<Me) = ^io g 

ip- (f?) = log 1 



¥) - 1 



^ cosh- 
?)-^cosh-(^ 



v 



(19) 



Thus the constants are determined by the initial data. Moreover, to satisfy the initial conditions 
for colliding plane waves, it can be shown that k± must satisfy 1 < \k ± \ < y2. It can then be 
seen that the initial wave in region II is 



*4+ 



2^2(1 - 



6{u) 



4(1 + 2 



^ 2 3V2a 2 k + 6ak 2 + t 2\/2A; + (l - k 2 s 



which, as behaves as 



(a + k. 



4+ 



(a + k 



^ (1-0 



f 2 



+ 



(1-^)3/2 



!i + £) s 



This clearly diverges unless a + k + has the values 0, +1 or —1. The opposing wave similarly 
diverges as r] — > 1 unless a + /c_ has the values 0, +1 or —1. Thus, within this family of solutions, 
initial waves which are bounded as £ — * — 1 or 77 ^ 1 and which give rise to solutions with a 
Killing-Cauchy horizon can only occur in a very limited number of possible cases. In fact, the 
only possibility up to a sign occurs when a = — 1 and k + = k_ = 1. This is precisely the 
degenerate Ferrari-Ibahez solution with a horizon |12j . 

Feinstein and Ibahez ^j] have also shown that the above result can be extended to solu- 
tions in which ip is given by (|18|) plus an arbitrary regular function that can be expressed as 
Jf^, a(u) Jo(uia)e luJ Pduj as in (|16j) . (This is currently the widest explicitly known class of col- 
liding plane wave solutions with a horizon.) It follows that perturbations which correspond to 
variations in the arbitrary function a{u) do not alter the character of the horizon at a = 0. 
On the other hand, perturbations of any of the parameters a or k ± lead to the occurrence of 
a curvature singularity in region IV and also give rise to to divergencies in the initial waves. 
Perturbations corresponding to additional components of b{uj) in (|1H|) similarly introduce both 
a curvature singularity in the interaction region and divergencies in the initial waves. 



6 The stability of solutions with horizons 

From the above discussion, it may be concluded that some perturbations which indicate that 
horizons in colliding plane wave space-times are unstable correspond to unbounded perturbations 



of the initial data. On the other hand, such horizons appear to be stable with respect to per- 
turbations which correspond to arbitrary variations in the function a(uj). In this case, however, 
although the corresponding perturbations of the initial waves are bounded, this corresponds to 
a situation in which the perturbation of one wave is exactly correlated with the perturbation of 
the other. Thus, it is not possible to perturb the initial data on one characteristic while keeping 
constant those on the other and at the same time preserving the character of the horizon in the 
interaction region. It follows that, for the above explicit solutions at least, the horizon must be 
unstable with respect to small independent perturbations of the initial data. 

To demonstrate the generic instability of any colliding plane wave solution with a horizon, 
consider a perturbation for which the initial data is given by 

MC) = + i>- (v) = ^b- (v) , (20) 

where ■)/>&+(£) an d ipb-iv) are the initial data which lead to a background solution with a horizon, 
and R+(£) is a perturbation function which is bounded for £ £ [— 1, 1] and such that R + £ — > 
as £ — ► —1. i.e. a perturbation has been included in one initial wave only and, when R+(£) = 0, 
a horizon occurs in the interaction region. 

Use can now be made of the solution of the characteristic initial value problem for colliding 
plane waves given by Hauser and Ernst [3] and implemented in [Jj. However, it is appropriate to 
present this in a slightly modified form which is more clearly related to that of the general non- 
linear case described in |Sj. In this method, it is first necessary to calculate the initial spectral 
amplitude functions A + (( + ) and from the initial data functions and ip-{rj) using 

the integrals 

A + (C + ) = j^— C ^ dg, A_(C-) = / f C ~ ?-« drj. (21) 

The solution in the interaction region can then be expressed as 

/' >™^±u «*)di + r ^^ i-(c.)d(.. (22) 

h v C+ - W(+ - v J—i -C-Vv- C- 

This method, which is based on the Abel transform, is linear. Thus, when the perturbation 
vanishes (i.e. when = 0), the resulting solution would be just the background solution 

i^bi^v) which contains a horizon. However, the perturbation term in (|2Ujl would lead to an 
additional perturbation component in the initial spectral data given by 

a p+ (c + ) = — f 1 -7=4^, v(c) = o, 

and this would lead to the solution in the interaction region being given by 



iP = ^ b + where V P = / V 4 / V(C+) ^C+- (23) 



- 1 yT^c:yc7+T 

It can then be shown that the scalar invariant + ?>^\ diverges as a — > (actually, it is 

the component ^2 which diverges), unless ip p ^ — ip PtTt — > in this limit. (This is equivalent to 
the condition (fTTjl .) 

However, the perturbation R+(£) is arbitrary and it is therefore possible to consider the com- 
ponent A p+ (C+) to be an arbitrary bounded function satisfying lim^ + _».i (a/1 — £ + A p+ (( + )) = 0. 
Moreover, it can be seen that 



(C+ - v) 



dC + , (24) 



which does not necessarily approach zero as a — > for an arbitrary choice of A p+ (( >+ ). (Indeed, 
it may well be unbounded, but this stronger result is not needed.) It can therefore be concluded 
that a perturbation of one initial wave will generically lead to a scalar polynomial curvature 
singularity in the interaction region in place of the horizon, i.e. the Killing-Cauchy horizon is 
unstable with respect to the perturbation (|2()|) . 



7 Conclusion 



Colliding plane wave solutions which contain a Killing-Cauchy horizon, only occur in very 
particular circumstances. Moreover, colliding plane wave space-times which are generated by 
bounded initial waves (apart from impulsive components on the shock fronts) and possess a 
Killing-Cauchy horizon form a highly restricted subclass of these solutions. It has been shown 
above that the only permitted perturbations of this subclass which preserve such a horizon are 
those in which perturbations of the two approaching waves are strongly correlated. However, 
when the initial waves are perturbed independently, then the horizon in the interaction region 
is necessarily replaced by a scalar polynomial curvature singularity. It is therefore concluded 
that the Killing-Cauchy horizon which sometimes occurs in colliding plane wave space-times is 
unstable with respect to general bounded perturbations of the initial waves. 

The above discussion has concentrated entirely on the "linear case" in which the approach- 
ing gravitational waves have constant and aligned polarization. However, it is expected that 
the behaviour of the nonlinear case in which the initial waves have variable or non-aligned po- 
larizations, or include electromagnetic waves, would be qualitatively similar. This expectation 
is supported by the arguments of Yurtsever fTI] for the vacuum case, but further analysis is 
required. 
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